function [y] = hmm_ti(a)


t1=a(1);
t2=a(2);
t3=a(3);
t4=a(4);
t5=a(5);
LPEM_TM_Lamellar_Ti;

y= sum(AA);

    y=-y/300;

savefig('figure ti_2')
end
